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The electronic properties of quarter-filled organic materials showing spin-Peierls tran- 
sition are investigated theoretically. By studying the one-dimensional extended Peierls- 
Hubbard model analytically as well as numerically, we find that there is a competition 
between two different spin-Peierls states due to the tetramized lattice distortion in the 
strongly correlated regime. One is accompanied by lattice dimerization which can be in- 
terpreted as a spontaneous Mott insulator, while the other shows the existence of charge 
order of Wigner crystal-type. Results of numerical density matrix renormalization group 
computations on sufficiently large system sizes show that the latter is stabilized in the 
ground state when both the on-site and the inter-site Coulomb interactions are large. 
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Quasi-one-dimensional (ID) organic conductors exhibit a variety of electronic states with spa- 
tially inhomogeneous charge, spin, and lattice structures. A typical example is the family of 
TMTSF 2 X and TMTTF2X, 1 ) where X denotes different anions. Recently, the existence of charge 
ordering (CO) has been identified in several of such quasi-lD compounds, which has renewed in- 
terest in these systems. It was first found in DI-DCNQ^Ag, a member of the -R1-R2-DCNQI2A 
family, with X= Ag or Li, and R\ and R2 taking different substitution groups to modify the DCNQI 
molecule itself, in which the charge pattern has been identified as a Wigner crystal-type one. 2 ' 3 -* 
Independently, analogous CO was predicted theoretically to also exist in TMTTF2A" based on the 
result of mean field calculations, 4 ) and soon after it was confirmed experimentally. 5,6 ^ 

In the two families menioned above, the ID 7r-band is quarter-filled in terms of electrons 
or holes. The wave vector along the chain direction for the CO state mentioned above is 4/cp, 
corresponding to the period of two molecules, which suggests the origin of this phenomenon to be 
the long range nature of Coulomb interaction as theoretically studied in the past. 7 ) Further studies 
have been performed on explicit models appropriate for the description of the electronic properties 
of DCNQI/TMTTF molecules, i.e., ID extended Hubbard models of quarter filling with on-site 
and nearest-neighbor Coulomb interactions, U and V, respectively. 4, 8 ~ n ) In those models, the CO 
insulating ground state is actually stable, in general, when both U and V are appreciably large 
compared to the transfer integrals. 
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Fig. 1. Schematic views of the ground states in quasi-one-dimensional organic conductors. Coexisting states of (a) 
charge ordering and Neel ordering, (b) charge ordering and spin-Peierls lattice distortion, and (c) dimer Mott 
insulating state and spin-Peierls lattice distortion. The grey area represents where the charge localizes, and 
the arrows and the ellipses represent the ordered spins and the spin singlet formation, respectively. The lattice 
distortions are also shown schematically by the thickness of the bonds. 

In the above compounds showing CO, below the CO transition temperature which is typically 
of the order of 100 K, magnetic phase transitions take place, to either antiferromagnetic (AF) or 
spin-Peierls (SP) state. These can be understood since in this CO state each charge is localized 
on every other site so that the spin degree of freedom acts as a quasi-lD spin 1/2 system. Then a 
competition arises between the AF state stabilized by the interchain exchange interaction and the 
SP state due to the ID instability coupled to the lattice degree of freedom. 12 ) In fact, the mean 
field calculations mentioned above, 4 ' expected to be relevant for the former case with sufficient 
three-dimensionality, show the coexistence of CO and AF for the ground state when V exceeds a 
critical value. This state is schematically shown in Fig. 1 (a), which is consistent with the low 
temperature AF states observed in DI-DCNQI 2 Ag, 2 ) and TMTTF 2 X with X = Br and SCN. 1 ' 13 ) 

On the other hand, in a recent experiment on TMTTF 2 X with X = PF 6 and AsF 6 , 6 ) the CO 
is actually found to also persist in the SP phase with lattice tetramization. Then this coexistence 
state can be schematically drawn as in Fig. 1 (b), which we call the CO-SP state. However, existing 
theoretical studies, mainly numerical ones, 14 ' 15 ) have failed to reproduce this CO-SP state within 
the so-called ID extended Peierls-Hubbard model, which is a natural extension of the extended 
Hubbard model to include the electron-lattice coupling. In these studies, only one other SP state is 
found, in which the lattice dimerization and the lattice tetramization both occur, as schematically 
shown in Fig. 1 (c). The lattice dimerization produces a spontaneous 'dimer Mott insulating' (DM) 
state, as will be discussed later, so we call this state the DM-SP state. 16 ) This state is relevant to 
several DCNQI compounds 2 ' 3 ) and to a classical example of SP compounds, MEM-(TCNQ) 2 . 17 ) 

In this paper, we will show that the coexistent state of CO and SP is stable in the strongly 
correlated regime, i.e., with large U and V, in contrast to previous studies. 14,15 ) This explains 
naturally the experimental facts noted above, and provides a unified view for the ground-state 
properties of the strongly correlated quarter-filled ID organic compounds on the basis of the CO 
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phenomenon. 

The ID extended Peierls-Hubbard model as in refs. 14 and 15 is expressed as 
H = 52t(l + VH)(cl* +lff + h.c) + fl>? 

i,cr i 

+U^2n i ^n il + V^2n i n i+1 , (1) 

i i 

where a is the spin index which takes f or j, and c- (cj CT ) denote the number operator and 
the creation (annihilation) operator for the electron of spin a at the ith site, respectively, and 
n-i = nj| + riii- U and V are the on-site and nearest-neighbor Coulomb energies, respectively. The 
electron-lattice interaction is included in the Peierls-type (or sometimes called the Su-Schrieffer- 
Heeger-type) coupling where the intermolecular motions result in the change of the transfer integral 
from t to t(l + m), while it expends the elastic energy Kuj/2, with K being the spring constant. 
Here Uj's are the normalized displacements which we treat as classical variables, and these are to 
be determined self-consistently to minimize the free energy of the system. 

First, let us analyze this model within the bosonization scheme 18, 19 - ) following the treat- 
ment of Yoshioka et al. 10 ^ The lattice distortion can be parametrized as u\ = Ud cos (ttx j /a) + 
ut cos (iTXi/2a + Xt), where X{ and a are the position of the ith site and the lattice con- 
stant, respectively, and xt is a phase factor which is determined in the following. We only 
consider the lattice dimerization Ud and lattice tetramization ut, which will be confirmed af- 
terwards in the numerical calculations that only these modulations are relevant in the ground 
state. The resulting phase Hamiltonian for the low-energy properties of this model is given by 
H = H P + H a + Ki/4 + H d + H t + H eh where 



(2) 



H° = £ J dx {j^ (<^+) 2 + K ° (W 2 } > (3) 
ft 1/4 = 9i/ a J dxcos46> + , (4) 
Hd = ~9d Ud J dxsin26> + , (5) 

ft* = ~9t u t J dx cos (e + -Xt-'j) cos <p + , (6) 

and H e i is the elastic energy term. 9 + and 4> + are the phase variables for the local density fluc- 
tuations of the charge and spin with a long wavelengh, respectively, and lie = —d x Q-/2ir and 
II^ = —d x (j)-/2ir are the corresponding "momenta" satisfying [9+(x),Ho(x')] = id(x — x') and 
[4> + (x),H ( p(x')] = i5(x — x'). Here, we neglected higher order nonlinear terms for the phase vari- 
ables as in ref. 10, since these terms have higher scaling dimensions. 
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The first two terms, 7i p and H a , describe the Tomonaga-Luttinger liquid where v p ,K p and 
v a , K a are the so-called Tomonaga-Luttinger parameters for the charge and spin parts, respectively. 
The SU{2) symmetry of the nondistorted model with m = requires K a = 1. 

The next two terms, H1/4 and Hd, are the quarter-filled and the half-filled Umklapp scattering 
terms, which are potential energies in favor of the phase variable to be fixed at 9 + = (or it/ 2) 
corresponding to the CO state and at 6 + = it/4 to the DM state, respectively. n ) The coupling 
constants are calculated in refs. 10 and 11 perturbatively from the weak-coupling regime as 5^/4 oc 
U 2 (U — 4V)/t 2 and oc U — AU(U — 2V)/t where A is some numerical constant, thus they 
become prominent in the strongly correlated regime. The DM state is analogous to that in strongly 
dimerized quarter-filled organic compounds such as K-ET2X, where the dimerization results in an 
effectively half-filled band. 20 ) However we should note that our DM state is possible only when the 
spontaneous dimerization Ud is finite, 21 ) thus it is a consequence of both the strong correlation and 
the ID instability toward lattice distortion. 

The Tit term is derived directly from the kinetic energy term in the presence of the lattice 
tetramization ut, where the coupling constant gt is proportional to i. 19 ) In the small U/t and V/t 
limits, this term produces the conventional weak coupling 2/cp-CDW state due to Peierls instability, 
which fixes the phase variables at 6 + = ir/4,<f> + = 0, and xt = 0. 14,19 ' 22 ) 

Therefore in the case of large U/t and V/t limits, which we are interested in, there occurs a 
competition between H1/4 and Hd, i.e., between the CO state and the DM state. Once the phase 
variable 9 + is fixed we may substitute its value into H, which result in an effective spin Hamiltonian, 

H' a = H a - g t ut J dx cos 4> + + H e i, (7) 

where we have optimized \t so as to gain the energy the most, i.e., \t = tt/4 for the CO state 
and xt = tt/2 for the DM state. This Hamiltonian is identical to the phase Hamiltonian for the 
SP problem in the ID Heisenberg model, where the spin singlet formation due to lattice distortion 
occurs even for infinitesimal spin-lattice coupling. 12 ) Thus in our case, the lattice tetramization 
ut always occurs once the electron-lattice coupling (and consequently, spin-lattice coupling) is 
included, then the phase variable of the spin is fixed at 4> + = and the spin gap opens. Namely, 
the CO state and the DM state result in the CO-SP state and the DM-SP state, respectively. We 
note that additional 2/cp-CDW appears in the presence of finite Ut where the order parameter is 
proportional to cos(2kpXi + 0+) cos</> + , 19 ) which is equal to cos(-7u/2) and cos(-7u/2 + tt/4) for the 
CO-SP state and the DM-SP state, respectively. 

To compare the relative stability of these states, the above bosonization procedure is not 
appropriate, not only since the coupling constants are obtained perturbatively from the weak 
coupling regime, but also because the treatment of the three nonlinear terms together with the 
lattice degree of freedom is a very subtle problem. Instead, we use the numerical density matrix 
renormalization group (DMRG) method 23 ) directly to model (1), which is essentially exact. The 
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periodic boundary condition is adopted and the lattice distortions Ui are calculated self-consistently 
as in ref. 14. We take the number of the states kept in the DMRG method, m, up to 250 and check 
that the m-dependence is very small in this region. 

The obtained results are summarized in Fig. 2, which shows the ground-state phase diagram 
on the plane of U/t and V/t for fixed values of 1/K = 1. For comparison, the phase diagram for the 
case of 1 / K = corresponding to the purely electronic extended Hubbard model 8 ) is also shown, 
where there are two phases, the Tomonaga-Luttinger liquid metallic phase and the CO insulating 
phase. 8 ' 10 ) The presence of the finite electron-lattice coupling 1/K = 1 makes three phases appear, 
the weak coupling CDW state, the DM-SP state, and the CO-SP state. One can see that the 
CO-SP state is stabilized in a rather wide range of parameters. We have noticed that the relative 
values of the ground-state energies for these states, and consequently their phase boundaries, are 
rather sensitive to the cluster size ./V when it is small such as ./V ^ 20 as in refs. 14 and 15. Thus we 
show in Fig. 2 the result for iV = 36 where such finite size effect is almost negligible, which provides 
qualitatively different results from those in the literature, especially in terms of the stability of the 
CO-SP state. 
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Fig. 2. Phase diagram of the extended Peierls-Hubbard model calculated by the DMRG method in the plane of U ft 
and V/t, for fixed values of 1/K = (dotted line) and 1 (filled line). For definitions of CO-SP, DM-SP, and CDW 
phases, see text. 

To observe the property of the CO-SP state in more detail, we show the lattice tetramization 
ut (ud is 0), and the electron density on each site in the CO-SP state as a function of V/t for a fixed 
value of U/t = 6 in Fig. 3 (a), and as a funcition of U/t for fixed V/t = 5 in Fig. 3 (b), both with 
1/K = 1. The electron density is parametrized as (rtj) = 1/2 + n 4 fc F cos(7TXj/a) + ri2k F cos(irxi/2a), 
where n^p and ri2k F are the order parameters for the 4/cp CDW, i.e., the Wigner crystal- type 
CO and the 2kp CDW, respectively. One can see that in the both cases of increasing the value 
of V/t and U/t, increases while ut and n 2 k F decrease. These results can be interpreted as 
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follows. It is natural that the degree of CO is enhanced when the degree of correlation, U/t or V/t 
is increased, thereby increasing n^k ¥ - In such a case, the decrease in u t and n2k F , representative 
of the spin singlet formation due to the SP state, can be explained by the strong U/t limit, where 
the effective spin exchange is proportional to the expectation value for both the nearest-neighbor 
sites to be occupied, (njnj+i). 24 ) Then the spin singlet formation energy decreases as the tendency 
toward CO is increased, since in the CO state the electron tends to occupy every other site. This 
result that and ut vary in an opposite way is consistent with the experimental observation 
that the CO and SP states 'compete' with each other, which is deduced from the variation of their 
transition temperatures as the pressure is applied. 6 ) 
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Fig. 3. Plot of the lattice tetramization ut, and the order parameter for CO and 2&f CDW, ri4k F and ri2fc F , respec- 
tively, for the CO-SP state, as a function of (a) V/t for a fixed value of 1/K — 1 and U/t = 6, and (b) U/t for a 
fixed value of 1/K = 1 and V/t — 5. 

The values of U/t for the actual compounds are believed to be approximately 5 ~ 7. 1 ' 2 ) There- 
fore the critical value of V/t ~ 4 in our calculations necessary to stabilize the CO-SP state is 
apparently very large. Actually, a reliable estimate yields V/t of about 2~3, 25 ) which indicates 
the region of DM-SP in our phase diagram. However, this estimated value for V/t is the Coulomb 
interaction between the neighboring molecules of the intrachain. The interchain Coulomb inter- 
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action should also be large, since the distance between the chains are rather close, especially in 
TMTTF compounds, 1 ^ although the interchain transfer integrals are small due to the anisotropy 
of the 7r-orbital. Thus, one should consider the parameter V/t in our ID model to be an effective 
one, which is possibly large enough to stabilize the CO-SP state. 

Finally, let us discuss some effects which we have not included in our model. Coupling to 
anions 14 ^ and/or the so-called electron-molecular vibration (e-mv) coup ling, 15 ) both of which result 
in the Holstein-type electron-lattice coupling modulating the on-site energy, help the stabilization 
of the CO-SP state as is studied in refs. 14 and 15. However, these effects should be small 
compared with the Coulomb interactions and also with the Peierls-type electron-lattice coupling 
due to the molecular displacements which we have studied. Therefore we believe that these effects 
are secondary. 

As for the alternation in the transfer integral which exists in TMTTF compounds, 1 ) as t(l + 
Ui) — > t(l + (— l)*Ad)(l + Ui), where is the degree of alternation, we expect it to be slightly 
disadvantageous for the CO-SP state. This is because the alternation results in modifying Hd in 
Eq. (5), as Ud — > Ud + Ad for the first order, and the gain in the energy by H compared to the loss 
in the lattice elastic energy S«n^ becomes larger for the DM-SP state. 

On the other hand, the value of 1/K has been chosen to be fairly large in the DMRG calcula- 
tions, in order to stabilize the numerical convergence. When we choose realistic values of 1/K, the 
phase boundary in Fig. 2 approaches that for 1/K = 0. This indicates the stability of the CO-SP 
state. However, this again is a subtle problem so numerical calculations should be pursued to be 
conclusive. 

In summary, we have argued that the coexistence of charge order and spin-Peierls lattice 
distortion observed in quasi-lD organic compounds is naturally reproduced by the effects of on- 
site as well as nearest-neighbor Coulomb interaction, together with the electron-lattice coupling 
of Peierls-type. Our results show that there is a competition between another spin-Peierls state 
coexisting with the spontaneous dimer Mott insulating state and that showing the existence of the 
charge order of Wigner crystal-type. The parameters for the actual compounds seems to be in the 
region near the phase boundary between the two states, so that by applying external fields such as 
pressure the ground state is expected to vary, from one to another, which is actually observed in 
TMTTF 2 AsF 6 . 6 ) 
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